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ABSTRACT 

In this contribution we discuss some of the main problems in high energy astro- 
physics, and the perspectives to solve them using different types of "messengers" : 
cosmic rays, photons and neutrinos. 

1. Introduction 

The birth of high energy astrophysics can be traced to nearly a century ago, when 
the balloon flights of Victor Hess established that a form of ionizing radiation, that 
was soon given the name of "cosmic rays", was arriving from outer space. Soon it was 
demonstrated that these "rays" are in fact ultrarelativistic charged particles, mostly 
protons and fully ionized nucleJB, with a spectrum that extends to extraordinary high 
energies. The existence of a large flux of ultrarelativistic particles arrived completely 
unexpected, and its origin remained a mistery that only now is beginning to be 
clarifled. The main reason for this very long delay in developing an understanding 
of this important physical phenomenon is that cosmic rays (CR) do not point to 
their production sites, because their trajectory is bent by the magnetic fields that 
permeates both interstellar and intergalactic space. 

Today we know that our universe contains several classes of astrophysical objects 
where non-thermal processes are capable to accelerate charged particles, both leptons 
(e''^) or hadrons (protons and nuclei), to very high energy. These relativistic particles, 
interacting inside or near their sources, can produce photons, and (in case of hadrons) 
neutrinos that then travel in straight lines allowing the imaging of the sources. At 
the highest energy also the magnetic bending of charged particles can become suffi- 
ciently small to allow source imaging. A detailed study of the "high energy universe" 
can therefore in principle be performed using three different "messengers" : photons, 
neutrinos and the cosmic ray themselves. This "multi-messenger" approach is still in 
its infancy, but has the potential to give us deep insights about the sites and physical 
mechanisms that produce these very high energy particles. 

Many of the proposed (or detected) acceleration sites are also associated with the 
violent acceleration of large macroscopic masses (one example is Gamma Ray Bursts), 
and therefore a fourth "messenger" : gravitational waves has the potential to give very 

^ Electrons have a steeper energy spectrum and constitute a fraction of only few percent of the 
flux; small quantities of positrons and anti-protons are also present. 
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important and complementary information about these astrophysical environments. 



2. The Cosmic Ray Spectrum 

The spectrum of hadronic CR observed at the Earth is shown in fig. [1], it spans 
approximately 11 order of magnitude in energy up to ~ 10^° eV. The CR spectrum 
is remarkably smooth and can be reasonably well approximated by a power law form 
(0 oc E~°') with a nearly constant slope. The most prominent spectral features are 
the "Knee" at ~ 3 x 10^^ eV, where the spectrum steepens from a slope a ~ 2.7 
to a ~ 3, and the "Ankle" at ^ 10^^ eV, where the spectrum flattens to have 
again the slope a ~ 2.7. At ~ 6 x 10^^ eV there is now clear evidence (from the 
HiRes^' and AugeJ^J collaborations) of a sharp steepening of the spectrum. Such a 
spectrum suppression had been predicted more that 40 years by Greisen, Zatsepin and 
KuzmirPl (GZK) as the consequence of the interactions of high energy particles with 
the photons of the 2.7 Kelvin Cosmic Microwave Background Radiation (CMBR). The 
non observation of this effect generated intense interest and an abundant literature. 
The observed spectrum suppression is consistent with the shape predicted for the GZK 
mechanism (pion production inp7 interactions), but other dynamical explanations are 
possible including fragmentation of heavy nuclei, or an intrinsic high energy cutoff 
in the CR sources. An additional (less evident) steepening of the spectrum (with 
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Figure 1: General structure of the CR energy spectrum plotted in the form (j>{E) E'^ versus logi? 



the slope increasing from a ~ 3 to approximately 3.3) the so called "2nd-Knee" at 
E 4 X 10^^ eV has recently emerged as potentially very significant. Explanations 
for the origin of the "Knee" , and the "2nd-Knee" and the "Ankle" are still disputed, 
and are considered as of central importance for an understanding of the CR. 

It is now knowiJ^ that most of the CR that we observe near the Earth are galac- 
tic", they are produced by sources inside our Galaxy, and are confined by galactic 
magnetic fields to form a "bubble" (or "halo") around the visible disk of the Milky 
Way with a form and dimension that are still the object some dispute. Extragalactic 
space should also be filled by a (much more tenuous) gas of CR, produced by the 
ensemble of all sources in the universe during their cosmological history. 

All galaxies (including ours) should be sources of CR, as the particles that "leak 
out" of magnetic confinement are injected into extragalactic space, however it is 
possible (and predicted) that there are additional sources of extragalactic CR that 
are not normal galaxies like our Milky Way, but special objects like for example 
Active Galactic Nuclei (AGN) or colliding galaxies. The extra-galactic component is 
expected to have a harder energy spectrum, and should emerge above the foreground 
of galactic CR at sufficinently high energy. The determination of a "transition energy" 
that marks the point where the extragalactic CR component becomes dominant is a 
very important problem (see discussion in section 4.3). 



3. Galactic Cosmic Rays 

Ultrarelativistic charged particles produced inside our Galaxy remain trapped by 
the galactic magnetic field, that has a typical strengtliB -B ~ 3 /iGauss. The Larmor 
radius of a charged particle in a magnetic field is: 

E 1.08 / E \ /iCauss 

rLarmor=^^-^ ^qIS ) " ^^pC (1) 

(one parsec is 3.084 x 10^*^ cm, or approximately 3 light years). Magnetic confinement 
becomes impossible when the gyroradius is comparable with the linear dimensions of 
the Galaxy. This corresponds to energy: 

E>ZeB i?Haio - 2.7 x 10^^ Z (:^^—] (^^] (2) 

\3 /iGauss J \10 Kpc J 

Most CR particles have much lower energy and a gyroradius that is much smaller than 
the galactic size. The motion of these particles can be well approximated as a diffusive 
process controlled by the random component of the galactic field. The time needed 



^ The magnetic field receives approximately equal contributions from a "regular" component 
(that in the galactic plane has field lines that run parallel to the spiral arms) and an irregular 
component generated by turbulent motions in the interstellar medium. Our knowledge of the field 
has still significant uncertainties and is very important fields of research. 
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for a CR particles to diffuse out of the galactic halo depends on its rigidity (indicated 
in the following for ultrarelativistic particles as E/Z). In zero order approximation 
(a spherical homogeneous halo) the confinement time of a CR can be estimated as: 

where -Rnaio is the halo radius and D{E/Z) is a rigidity dependent diffusion coefficient. 
An important consequence is that (assuming stationarity) the population Na{E) of 
CR of type A and energy E in the Galaxy is: 

Na{E) - Qa{E) r{E/Z) (4) 

where Qa{E) is the injection rate and t{E/Z) is the (rigidity dependent) residence 
time. Therefore the observed CR spectrum has not the same energy distribution of 
the particles near their sources, but is distorted and steepened by confinement effects: 

The determination of the shape of the "source spectrum" Q{E) and of the length 
and rigidity dependence of the confinement time are obviously important problems. 
A powerful method is the measurement of the spectra of "secondary nuclei" (such 
as Lithium, Beryllium and Boron) in the CR. These nuclei are very rare in normal 
matter, but they are relatively abundant in the CR flux, because they are produced 
in the spallation of parent nuclei such as Carbon and Oxygen, a process where the 
nuclear fragments maintain the same velocity and therefore the same energy per 
nucleon as the incident nucleus. The energy spectrum of secondary nuclei is expected 
to be proportional to the square of the confinement time (its own, and of the parent 
particle). The measurement J^J indicate that while "primary" nuclei have spectra 
oc E""^-"^ , secondary nuclei have steeper spectra oc E~^-^ . The conclusion is that the 
confinement time decreases with rigidity as r oc E^^ ~ E~^'^, and therefore the 
injection spectrum of CR at their sources has a slope ainj — Oobs — ^ — 2.1. The 
absolute value of the confinement time can be estimated from the relative amounts of 
secondary and primary nuclei in the CR flux, with additional information obtained 
studying the abundance of unstable nuclear isotopeJ^ with appropriate lifetime (in 
particular the nucleus Beryllium-10 with a half like of 1.5 Million years) and is of 
order of 10 million years for a rigidity of a few GeV. 

The power law rigidity dependence (r oc E~^) of the confinement time, and of the 
diffusion coefficient {D oc E^) is not unexpected, because for a turbulent magnetic 
field one predicts a power law form for the diffusion coefficient with a slope that 
is a function of the power spectrum of the field irregularities. For the theoretically 
favoured Kolmogorov turbulence one expects t{E) oc E~^/^. This results would imply 
a steeper injection spectrum oc E~'^'^'^ . The discrepancy between this theoretical 
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prediction and data is is still under discussion. 



3.1. The "SuperNova Paradigm" for galactic Cosmic Rays 

A confinement time of a few million years is in fact very short with respect to the 
age of the Galaxy, and therefore CR must continuously be produced in our Galaxy 
(that is in a stationary state with approximately equal CR injection and loss). The 
power of the ensemble of CR sources in the Milky Way can be estimated as the ratio 
between the total energy of CR in the Galaxy divided by their average confinement 
time: 

^Milky Way ^ Pcr{x&) Vcff ^ 2 X 10"^^ ( — \ (5) 
(Tct) \ S / 

For this numerical estimate we have used a local CR energy densitjE Pct{xq) = 
1.6 eV/cm^, an effective volume of 170 Kpc^ and a average confinement time of 
20 Million years. This surprisingly large power requirement (approximately 50 mil- 
lion solar luminosities), is a very important constraint for the identification of the 
galactic CR sources. In fact, s imp ly on the basis of this energy condition SuperNova 
(SN) explosions were proposed-^ as the most likely CR source. The average kinetic 
energy of a SN ejecta is of order L^in — 1.6 x 10^^ erg (this corresponds to the ki- 
netic energy of 10 solar masses traveling at f ~ 4000 Km/s) . For a SN rate in our 
Galaxy of order every (30 years)"^ the power associated to all SN ejections is then 
L|™j^ ~ 1.7 X 10^^ erg/s, and it is therefore possible to regenerate the Milky Way 
CR if approximately 10-20% of the kinetic energy of the ejecta is converted into a 
population of relativistic particles. 

In addition to these energy balance considerations, in the 1970's a dynamical 
argument emerged in favour of the SuperNova hypothesis, when it was understood 
(see Drurji^ for a review) that the spherical shock waves produced in the interstellar 
medium by the (supersonically moving) SN ejecta can provide the environment where 
particles, extracted from the tail of a thermal distribution, are accelerated up to very 
high energy generating a power law spectrum with a well defined slope of order 
a ~ 2 + e with e a small positive number. That is precisely the injection spectrum 
needed to generate the observed CR. 

The basic concept behind this theory is an extension of the ideas developed by 
Enrico Fermff^, who in 1949 made the hypothesis that the acceleration of cosmic 
rays is a stochastic process, where each CR particle acquires its high energy in under- 
going many collision with moving plasma clouds. The clouds carry (in their own rest 
frame) turbulent magnetic fields and act as 'magnetic mirrors" imparting on average 



The estimate of the CR energy density is of the same order of magnitude of the energy density 
of the galactic magnetic field (for i? ~ 3 /iGauss one has pB = B"^ /{Sir) ~ 0.22 eV/cm'^). This is not 
a coincidence, and has a transparent dynamical explanation. The galactic magnetic field is mostly 
generated by electric currents transported by the CR, while the CR themselves are confined by B, 
and therefore the two quantities are close to equilibrium. 
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a positive "kick" to the scattering particle with {AE)/E oc (with jS the cloud 
velocity). It is easy to see that this process generates a power law spectrum, with 
higher energy particles having performed a larger number of collisions. 

The crucial new element introduced in the 1970's is the presence of the shock wave. 
Charged particles are now accelerated scattering against magnetic (in the plasma rest 
frame) irregularities present both upstream and downstream of the shock front that 
act as Fermi's clouds, or converging magnetic mirrors. The new geometry allows a 
more efficient acceleration [((Aii^)/i?)crossing oc /3shock] and the (mass, momentum and 
energy conservation) constraints of the fluid properties across the shock determine 
the slope (a ~ 2 + e) of the accelerated particle spectrum. 

This mechanism (diffusive acceleration in the presence of shock waves) can operate 
every time that one has a shock wave in an astrophysical fluid. The spherical blast 
waves of SN ejecta are one example of this situation but several other are known to 
exist. In fact shocks are generated every time that macroscopic amounts of matter 
move at supersonic speed. Particularly interesting case are the jets emitted by Gamma 
Ray Bursts, by accreting Black Holes of stellar mass (microQuasars) or supermassive 
(Active Galactic Nuclei). In all these objects there in in fact evidence for particle 
acceleration. 

The "SuperNova Paradigm" has a simple and very important implication: in 
the vicinity of young SuperNova Remnant (SNR) one should flnd a population of 
relativistic hadrons (protons an nuclei) with a spectrum close to the injection one 
(E-(2+^)) and a total energy of order ~ 0.2 x 10^^ erg. These relativistic particles can 
interact with the interstellar medium around the SN producing neutral and charged 
pions that then decay generating photons and neutrinos with a spectrum that to 
a good approximation follows the same power law of the parent proton spectrum. 
Describing the relativistic proton population as Np{E) ^ Kp E''^, and approximating 
the conflnement volume as homogeneous and with density Ut, the rate of emission of 
photons and neutrinos is approximately: 

N^(^^){E) = Kp {cappTiT) Zp^^(^){a) E~" (6) 

where app is the pp interaction cross sectioii^ and Zp^-y(^iy) are adimensional propor- 
tionality factors. For a proton spectrum oc that has equal energy per decade 
of energy the Zp^^i^^^ factors are well approximated by the (approximatgely energy 
independent) fraction of the energy of the interacting protons that goes into photons 
(or neutrinos), with 

^P'y ^'p^ 0.X5. Usm^ ttiG th-GorGtiCcLlly GxpcctGci slope oi — 2, 
flxing the constant Kp by the total energy contained in the relativistic particles, as- 
suming a spectrum extending in the interval [-Emin? -^max], a source at a distance d 

^¥or simplicity we use a notation that indicates only the dominant components for both the CR 
and the target medium. 
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and integrating above Eth the photon (or neutrino) integrated flux at the Earth is: 



. Elf Z, 



/Kpc\2 ( Ut \ ( El°' \ /TeV\ _2 „i 




0.9 X 10-1^ Kt^ cm-^ s-i (7) 

\ d J \cm--y I 1050 erg/ V ^th / 



whith (i the source distance. Equation ([7]) is a key prediction of the "SNR paradigm" 
for the acceleration of the bulk of cosmic rays that can be tested with photon obser- 
vations in the GeV/TeV regions, and soon also with neutrino telescopes. 

In 20o4^ the HESS TeV Cherenkov telescope observed the SuperNova Remnant 
RX J1713.7-394ffl as a very bright source of TeV photons. The property of the 
photons from this source, taking into account our knowledge of the density of the 
environment around the object are consistent!!^ with the expectation of the "SNR 
paradigm" for galactic cosmic rays. 

Several other young SNR have also been detected with TeV photons (this in- 
cludes the detection of Cassiopea , while for other objects one has only upper 
limits. These observation of SNR in high energy photons gives support to the "SNR 
paradigm" for CR acceleration (for a more critical view of the situation se^^). 
This conclusion is however not unambiguous; also for the best candidate source 
RX Jl 713. 7-3946 a "leptonic origin" of the radiation (inverse Compton of relativistic 
e"^ on the radiation fields around the SN) cannot be entirely excluded. 

It is remarkable that the neutrino flux that should accompany the photon flux 
in case of a "hadronic" origin (vr" — 77 decay, accompanied by vr"*" — >• fi'^ i'^ 
(e"*" z/e'F^) z/ + /i chain decay) is detectable, at the level of few neutrino induced up- 
going muon events per year, in a Km'^ size neutrino telescope placed in the northern 
hemisphere. Photon observations in the GeV region with detectors on satellites (Ag- 
ile and GLAST) have also the potential to test the "SuperNova paradigm" for CR 
production. 



3.2. The "Knee", and the maximum acceleration energy 

The mechanism of diffusive particle acceleration generates a power law spectrum 
only up to a flnite maximum energy -Emax- This maximum energy is determined by 
the product of the acceleration rate times the total time available for the acceleration. 
Simple considerations allow to estimate -Emax — Z eB R /3shock where B is the typical 
strength of the magnetic fleld, R the linear dimension of the acceleration site and 
Ahock the shock velocity. Note that this estimate is more stringent than a simple 

° An object discovered in X-rays by the ROSAT instrument and then associated, with a very 
high degree of confidence, to a supernova observed and recorded by Chinese astronomers in the year 
393 of the current era. 
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magnetic containement condition by a factor /^shock < 1- Substituting values that are 
relevant for a SNR one finds: 



E^,.^Ze5i?Ahock^0.4xlOi^eVZ f-^") ^ ] (^] (8) 

\o pc J y3/iGaussy y 0.03 J 

This estimated maximum energy is much smaller than the highest observed energies, 
and therefore it is clear than other type(s) of CR sources must exist. 

The estimate ([8]) is however not too dissimilar from the "Knee" energy (at E ~ 
3 X 10^^ eV, and it has been suggested that the "Knee structure" marks in fact the 
maximum energy for proton acceleration in SNR. The apparent smooth behaviour 
of the spectrum above the knee could hide a succession of cutoffs with increasing 
energies (-Emax = ^ -^max) heavier nuclear species. In this model the SNR spec- 
trum would extend up to the maximum energy for iron with E^^^ ~ 0.8 x 10^'' eV. 
The key prediction here is the gradual change (with increasing (A)) CR compositon 
above the knee, or more precisely a set of different "knees" at identical rigidity, and 
therefore scaling in energy proportionally to the electric charge Z. Unfortunately the 
measurement of particle mass in Extensive Air Showers in the region of the kneqj are 
not easy and suffer from significant uncertainties in the modeling of the shower (see 
below in section 6). The data of KascadJ22l ^^d other experiment are (with large 
errors) roughly consistent with this idea. 

Equation ([8]) for the maximum energy is valid in general, for all astrophysical 
environments where shocks are present in a diffusive medium, and setting /?shock = 1 
the equation simply gives a very general condition for containement in the source. 
Since one observes particles with energy up to 10^° eV, the acceleration of these 
particles requires sources with a sufficiently large B R product, that is sources with 
a sufficiently strong magnetic field and sufficiently extended to satisfy condition (jH]). 
This point has been stressed by Michael Hillas^ that initiated a systematic study 
of all possible astrophysical environments that are possible candidates for ultra high 
energy acceleration on the basis of equation ([8]). The bottom line of this analysis is 
that this very simple constraint is very "selective" and only few objects are viable 
for acceleration up to ~ 10^° eV. The main candidates are Active Galactic Nuclei, 
and the jet-like emission associated with the Gamma Ray Bursts. 

3.3. Alternative Models 

A possible alternative to the SNR paradigm is the "cannon ball model" of Dar and 
De Rujulsff^. In this model the Energy source for the production of CR is ultimately 
the same as in the SNR paradigm, and is the gravitational binding energy released 
during the collapse of massive stars at the end of their evolution, but the dynamics 

^ The best method is the simultanous measurement of the electromagnetic and muon component 
of the shower at ground level. 
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of CR acceleration is deeply different. In the cannon ball model the explosion due to 
a gravitational collapse produce in most of perhaps all of the cases not only a quasi- 
spherical mass ejection with an initial (non-relativistic) velocity of order 10^ Km/s, 
but also the emission along the rotation axis of the collapsing star of "blobs" of 
ordinary matter each having a mass of order 10^^ grams and an extraordinary Lorentz 
factor of order F = 1000; the "cannon balls" emission, according to the authors is 
responsible for most of the Gamma Ray Bursti^, with a mechanism that differs 
from the often discussed "Fireball model"^ '^ ^ ^^ ^ Cosmic rays are accelerated by 
the collisions of the cannon balls while they decelerates slowly colliding with the 
interstellar medium. 

In this work we do not have the space to discuss critically this model, see referencJ^H 
for criticism (and the authors' respons^^). 

4. Extragalactic Cosmic Rays 

Most ultrarelativistic particles injected in extragalactic space have negligible in- 
teraction probability, and only lose energy because of the universal redshift effects. 
The important exceptions are those particles that have a sufficiently high energy to 
interact with the abundant photons of the Cosmic Microwave Backround Radiation 
(CMBR). For protons one has to consider two important thresholds associated to pair 
production (p + 7 — > p + e~^e~): 

*^ 10 ^ ^ 

and pion production (p + 7 — >■ + vr): 

_ UhL^ ^ 6 X 10^9 eV (10) 

For heavy nuclei one has to consider the photodisintegration reactions with threshold 
(related to the process A + •y ^ {A — 1) + N): 

*^ 10 V8MeV/ ^ ' 

where Eb is the binding energy of the displaced nucleon. Note that the threshold 
energies for protons and iron nuclei {A ~ 56) differ by a factor 2-3. Particles above the 
thresholds (Q and (ITOl) lose rapidly energy cannot propagate only for a cosmologically 
short time. This results in a strong flux suppression. 

The (space) average extragalactic density {na{E)) of particle type a (for example 
protons, or neutrinos) can be calculated integrating the source emission over the 
entire history of the universe, taking also into account the effects of energy loss or 
absorption. A general expression for the average particle density is: 

roc p 

{n{E)) = j^ dt J dEgq{Eg,t) P{E,Eg,t) (12) 
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where the double integral is over the emission time and the emission energy, q{E, t) 
is the (space averaged) emissivity for the particle considered (that is the number of 
particles emitted per unit of time, comoving volume and energy at universal time t), 
and P{E,Eg,t) is the probability that a particle generated with energy Eg at time t 
survives until the present epoch with energy E. 

Expression (fT2l) can be simplified if one makes the assumption that the energy 
loss is a continuous deterministic process, and the probability P can be approximated 
with the expression: 

P{E,Eg,t)=6[E-Eg{E,t)] (13) 

where the function Eg{E,t) is the energy of a particle observed now with energy 
E after propagation backward in time to time t. This approximation is essentially 
perfect for neutrinos, and also sufficiently accurate for protons. 

One can rewrite equation (fT2ll changing the integration variable from t to the 
redshift z with the Jacobian factor: 

^ = I (14) 

where H{z) the Hubble constant at the epoch redshift z: 



H{z) = Ho ,JQ,r. (1 + zy + Qa + {1-^a- ^m) (1 + z)^ , (15) 

and performing the integration over Eg using the delta function. The result can be 
written in the form: 

(niE)) = ^ q{E,0) CiE) (16) 
where C(-^) is an adimensional "shape factor": 

The expression flTB]) has a trasparent physical interpretation: q{E, 0) tells us how 
many particles of energy E are injected per unit time and unit volume in the present 
universe, multiplying by the Hubble time Hq^ we obtain a first order estimate of 
the present density, and the "shape factor" contains all complications related to the 
expansion of the universe, the cosmological evolution of the sources, and the energy 
losses of the particle considered. 



4.1. The proton extragalactic spectrum 

We can now use expression flTB]) to compute the extragalatic proton flux assuming 
that it is related to the space averaged density by the simple relatior@ 

ME) = ^ ME)) = ^ ^ q{E, 0) CiE) (18) 

Sin general the effects of the magnetic field, and of the "granularity of the sources can result in 
modification of this assumption. 
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The simplest assumption for the injection is a power law form: q{E) = qq E~" 
that depends only on two parameters, the slope a and the absolute normalization. 
The calculation of the flux has three obvious steps, the first is the calculation of the 
evolution of the proton energy with redshift; the second is the calculation of the shape 
factor C{E) according to equation (fTTI) . and finally one can use equation (ITSl) . 

The results of a numerical calculation of Eg{E, z) are shown in fig. [2l The differ- 
ent curves correspond to the "backward in time" evolution of the energy of protons 
detected "now" with E between 10^^ eV and 10^^ eV. For each energy E one finds (at 
a certain redshift z) a sort of "photon wall", where the E{z) starts to grow faster than 
an exponential. Clearly at this critical redshift the universe (filled with a hotter and 
denser radiation field) has become opaque to the propagation of the protons. Note 



> 




Redshift (z) 

Figure 2: Redshift (time) dependence of the energy of protons observed with different final energy 
at the present epoch. 



that the time evolution of the energy of a proton is completely independent from the 
structure and intensity of the magnetic fields, because a magnetic field can only bend 
the trajectory of a particle and the target radiation field is isotropicEI. 

The result of the calculation of the "shape factor" C(-^) for ^ power law injection 
is shown in fig. [3l It is important to note that the function C{E) has a non trivial 
energy dependence, and while the injection spectrum is a smooth power law will show 



''Clearly in this case the redshift mantains its one to one correspondence with the time of emission 
of the particle, but loses its simple relation with the space distance of the particle source. 
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Figure 3: Shape factor C{E) for protons, calculated for a power injection with slope a = 2.0, 2.2, 2.4 
and 2.6. The solid lines are calculated for a constant value of the source power density, while the 
dashed lines assume for Cp{z) the redshift dependence shown estimated by Ueda^^ll ^j^g X-ray 
luminosity of AGNs. 



some structure, the "imprints" of the energy loss processes (pair and pion production) 
of the protons. For energy E ^ E^h^ , only adiabatic (redshift) energy losses are 
significant, and in this case the the shape factor ({E) tends to a constant value: 

c^fd^^d+.r (19) 

Jo J^{z) Lq 

{C{z)/Co is the cosmological evolution of the source luminosity). Accordingly, the 
power law injection results in an observable flux that is a power law of the same slope. 
At higher energy the energy losses due to the pair production and pion production 
thresholds (equations (Q and (fTOl) ) leave their "marks" on the shape of the energy 
spectrum. 

The results of the calculations for the obervable flux are shown in fig. H] and 
fig. [5l The calculations were performed with different slopes of the injection power 
law a = 2.0, 2.2, 2.4 and 2.6, and with two different assumptions for the evolution 
of the cosmic sources, a constant injection, and evolution equal to the one estimated 
ijpH fQj; AGN hard X-ray emission. The normalization of the calculation is 
chosen in order to have a good fit to the observed data in the very high energy region. 
Some comments about the results shown in fig. H] are in order: 
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E (eV) 

Figure 4: Extragalactic proton spectrum calcul ated for a power law injection with different slopes 
a, for no-evolution of and with the evolution oP^l, xhe data points are from the Tibet arrajEH 
and the HiRes detector^. 



(i) The shape of the observed flux above E > 5 x 10^^ eV is consistent with the 
presence of the suppression predicted by The implication of this point are 
discussed below in section 4.2. 

(ii) The shape of the calculated spectrum is determined by the slope a of the 
injection spectrum and also by the assumptions about the redshift dependence of the 
injection. Actually this redshift dependence is of small importance for E > 5x 10^^ eV 
because the lifetime of a very high energy particle is brief and they do not probe 
deeply into the past history of the universe. For lower energy the inclusion of source 
evolution is significant. Lower energy particles can arrive from a more distant past 
when the emissivity was higher, and therefore, for the same shape of the injection, 
the inclusion of source evolution results in a softer spectrum. As a rule of thumb, 
including for the injection the redshift dependence of the AGN luminosity (or also 
of the stellar formation rate) has an effect similar to the softening of the injection 
spectrum of order Aa ~ 0.2. 

(iii) It is remarkable that if one choses an injection spectrum with a shape oc E'"^'^ 
and no source evolution (or a source oc with cosmic evolution) with appropriate 
normalization, the calculated flux can provide a good match to the data in the entire 
region from E ~ 10^^ eV up to the highest energies. The possible implications of this 
result are discussed below in section 4.3. 
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Figure 5: High energy CR spectrum. For clarity only the data of the Tibet array^ (below E ~ 
IQi^ eV), HiReiH and Auger^* (at higher energy) are shown. The spectrum is plotted in the form 
(t>{E) versus logi?. On can note the different energy scales for the Hires and Auger experiment. 
The Tibet data points are calculated using three different interaction models for the reconstruction 
of the shower energy. The lines describe the shape of an extragalactic components with injection 
spectrum E~°' and normalized to fit the high energy data. 



4-.2. The "GZK controversy" 

For several years the "dominant question" in CR physics has been the possible 
existence of a large flux of particles with energy larger that 10^° eV. To solve this 
very surprising result two main arguments have been proposed. One possibility is 
the existence of an additional source of particles of extremely high energy from the 
decay of super-heavy (mass of order of the GUT mass Mqut — 10^^ eV). A second 
possibility is the existence of violation of the Lorentz invariance. The process jrj 
Nn that is at the basis of the expected suppression of the proton flux is in fact studied 
in the laboratory mostly in the p rest frame, while for CR it is relevant in a frame 
where the p has a Lorentz factor Ep/rUp ~ 10^^. One can therefore speculate that 
"small" violations of the Lorentz invariance could suppress the cross section for high 
energy CR. Ultra High Energy Cosmic Rays therefore offered the exciting possibility 
of being the key observational field for the exploration for new physics. 

The new results of the HiRei^^ and Augei^ collaborations demonstrating the 
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existence of a suppression of the CR flux at ii^ ~ 6 x 10^^ eV have however very sig- 
nificantly changed the "scientific landscape" around the field. The phenomenological 
motivations for "exotic physics" in UHECR have essentially disappeared. Of course, 
the theoretical arguments behind the existence of supermassive topological defects as 
relics from the early universe, and (perhaps more controversially) of possible viola- 
tions of the Lorentz invariance remain in existence, and therefore UHECR can be used 
to set limits, or look for hints of these phenomena. However the field of UHECR has 
now returned to its status as an "astrophysical problem" , where the main questions 
is to establish what are the sources of these ultrarelativistic particles, and what are 
the physical processes responsible for the acceleratioiiH. 

Does this new situation makes the study of CR less interesting? In a certain 
sense the answer is obviously yes. The significance of the "exotic physics" for the 
UHECR would clearly be of truly revolutionary significance. On the other hand the 
astrophysics of CR, or more in general study of the origin of the different forms of 
high energy radiation remain a field of fundamental importance. 

The development of fundamental Science has historically always been intimately 
connected with the discovery and understanding of astrophysical objects. Steps in 
this remarkable history have been the understanding of: the structure and dynamics 
of the solar system; the source of energy in the Sun and the stars; the quantum 
mechanical effects that sustain White Dwarf stars; the formation of neutron stars 
and their connection with SuperNova explosions. 

The sources of the high energy radiation in the universe (some of which are still 
undiscovered) are very likely (together with the early universe) the most "extreme" 
environments where we can perform the most stringent tests for physical laws in 
conditions that are often simply not achieavable in Earth based laboratories. In 
particular, the regions of space near compact objects or near the horizon of Black 
Holes are very likely particle acceleration sites. It is also possible (or perhaps likely) 
that the still unexplained cosmological Dark Matter will play an important dynamical 
role in precisely these acceleration sites (for example near the center of galaxies), and 
therefore one could have to disentangle the physics of particle acceleration and Dark 
Matter annihilation. In a nutshell: High Energy Astrophysics remains a crucial field 
for fundamental Science. 



4.3. Galactic to Extragalactic transition. 

Berezinsky and collaboratorJ^Hl have been the first to note the remarkable fact 
that the injection of a smooth (power law) spectrum of protons (with a suflicienly 

' A posteriori it could be an interesting problem to consider why essentially the entire community 
of CR physicists has remained so "captured" with the idea of the existence of "super-GZK" particles" 
on the basis of an evidence of weak statistical significance and with the possibility of significant 
systematic errors. The "temptation" of the possible discovery of a result of extraordinary importance 
probably played a significant role for both observers and theorists. 
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small contributions of helium and other nuclei) can describe the UHCR not only above 
the "ankle" but in the entire energy range E > 10^^ eV. These authors stress the 
point that the prediction depends on only few parameters (the slope a, the injection 
absolute normalization and the z dependence of its evolution) can provide a very 
good fit to the data, and have concluded that this is very unlikely to be a simple 
coincidence, and is evidence of a real effect. 

The "Ankle" feature in the CR spectrum has been for a long time assumed to mark 
the transition energy where the softer galactic component is overtaken by the harder 
extragalactic one. In the Berezinsky et al. model the flattening of the spectrum 
emerges from an entirely different physical mechanism. In the entire energy range, 
both above and below the "Ankle", the spectrum is due to a single, proton dominated 
component, injected with a smooth (power law) spectrum. The visible structure in the 
spectrum is interpreted as the imprint of energy losses due to e~*"e~ pairs production. 
This effects is mostly effective at an energy range around E ~ 10^^ eV, and form a 
"dip" in the spectrum. The energy range of the effect is controled by the kinematical 
threshold for pair production at low energy, and by the fact that the energy lost 
for the production of an e~^e~ pair is oc ml/ E decreases with E. The shape of the 
sprectrum distortion or "Dip" is in fact determined by simple, fundamental physics 
cosiderations and robustly predicted. The form of the spectrum is of course strongly 
dependent on the value of the slope a. When the slope is a ~ 2.6 (for no-evolution 
of the injection) or a ~ 2.4 (for AGN/SFR source evolution) the spectral shape can 
describe the data in the entire energy range: [10^^, 10^''] eV, 

The name "Dip model" encodes this new dynamical explanation and renames the 
structure traditionally called "Ankle" . In this model the transition between galactic 
and extra-galactic cosmic rays happens at lower energy and corresponds to the less 
prominent spectral steepening feature called "2nd-Knee" . For more discussion about 
observations of this structure seJ^SJ 

The "Dip Model" has attracted considerable interest, and is a leading contender 
for the description of the highest energy cosmic rays. In my view, the questions of 
where is the transition between galactic and extra-galactic cosmic rays, and if the 
"Ankle" is really a "Dip" or viceversa remain open, at the center of an important 
controversy, that will be soon settled by improved observations. The implications of 
this problem are important and wide. 

If the "Dip Model" is correct, life for the observers would be significantly simpli- 
fied. The reason is that it allows "calibrations" for the measurements of both energy 
and mass of the CR starting from E ~ 10^^ eV. It is clearly possible to use an ob- 
served feature of the spectrum (with a well understood physical origin) to determine 
the scale of the energy measurement^. In fact, the idea to use the "GZK suppres- 

j Matching the spectral shapes of different experiments one can estimate the ratios between their 
respective energy scales. To establish the absolute energy scale one obviously need to compare the 
observed spectra with a "template" based on an understood physical mechanism. 
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sion" to fix the CR energy scale is now 40 years old. In the "DIP model" one has 
a spectral feature (the "Dip") that allows to play the same game at lower energy 
and with much higher statistics. Berezinki and collaborators have already suggested 
corrections factors for the energy scale of the different UHECR experiment. It should 
be streessed that the vahdity of the "Dip Model" depends on the fact that the energy 
scale of the HiReS experiment is approximately correct, while the energy scale of the 
Auger experiment is significantly underestimated, as shown in fig. [5l 

The same discussion can be made for the CR composition measurements. An 
established "GZK" feature implies a proton rich spectrum. In the "Dip Model" most 
particles above E ~ 10^^ eV are protons. With this knowledge one can use the obser- 
vations on the shower development to extract information on hadronic interactions 
at high energy (see discussion in section 7). 

The relatively "soft" injection {a = 2.4-2.6, depending on the assumptions about 
the source resdshift evolution) of the "Dip model" is not expected in several accelera- 
tion models for both AGN's and GRB's that predict a slope closer to a = 2. It should 
be noted however that the injection considered here is the space averaged one, and it 
is possibl^^^ that the spectra of individual sources are all flatter (for example with 
a ^ 2 but have different high energy cutoffs, that combine to an average spectrum 
that can be approximated with an effective softer power law. This explanation is 
however problematic, and it is unclear if it can naturally result in spectrum with the 
required smoothness. 

The determination of the galactic/extra-galactic transition energy has obviously 
important consequences also for the properties of the galactic spectrum, and the 
galactic sources. In the models that have the transition at the "Ankle" it is necessary 
to assume that our own Galaxy contains sources capable to accelerate particles up to 
10^^ eV and above. In the "Dip model" with a lower energy transition one has a less 
stringent (by a factor 10 of more) requirement for the maximum energy of the Milky 
Way accelerators. 

The transition energy between galactic and extra-galactic cosmic rays obviously 
marks the passage from a softer to a harder component. It is therefore natural to 
expect that the transition appears as an hardening of the spectrum. The 'Ankle" is 
the only observed hardening of the spectrum, in the entire CR energy range, this has 
lead to its identification with the transition energy. In the Dip model the galactic 
to extra-galactic transition corresponds however to the "2nd-Knee", that is to a 
spectral steepening. This may appear, and indeed it is surprising. The situation 
where the transition to a harder component appears as a softening is in fact possible 
but is requires the existence of two condition, that while perfectly possible are, in the 
absence of a dynamical explanation surprising coincidences: 

1. Both galactic and extra-galactic components must have a steepening feature at 
approximately the same energy E*. The extragalactic component must however 
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be the harder of the two, both below and above E*. 
2. At the energy E* the two components must have approximately equal intensity 

(0gal(^*) ~ 0extra(^*)). 

Let us consider a situation where the galactic component has slopes Ogai and /3gai > 
ttgai below and above E*, while the extragalactic component has slopes Oextra and 
Pextr > ttextra- If the relative normalization of the two component is "about right", 
and if the 4 slopes are ordered as: 

Q^extra ^ Clgal ^ /^extra ^ Z-^gal 

then one will observe a spectrum that steepens from the slope ctgai to the slope Pextra., 
with the steepening marking the transition to the harder extra-galactic component. 
The two slopes Oextra (extragalactic component below the transition) and /5gai (galactic 
component above the transition) are not easily measured because they are "hidden" 
by the dominant component. 

The discussion above is just a toy model, but in fact it is also an approxima- 
tion of what happens in the "Dip model" . At the transition energy (coincident with 
"2nd-Knee" at -E* ~ 4 x 10^'' eV) one observes a spectrum softening from a slope 
a ~ 3.0 to a slope ~ 3.3. The spectrum below the knee reflects the shape of a more 
abundant galactic component (0gai oc E~^), while the harder (0extra E'"^-^) extra- 
galactic proton component is hidden below. Near the transition energy the galactic 
component two things happen: the galactic component steepens sharply (presumably 
because of the "end" of the galactic accelerators), while also the extragalactic com- 
ponent steepens because of the opening up of the pair production channel for energy 
loss, taking (just below the "2nd-Knee") the approximate slope ~ 3.3. The relative 
normalization of the two components is such that one observes a single steepening 
transition (approximated as the slope transition 3 — >■ 3.3). 

This situation is possible, but it is not natural, and seems to require some "fine 
tuning". In fact the same critical observations have already been formulated in the 
past discussing the "knee" (-Eknee ~ 3 x 10^^ eV) as a possible transition between two 
different CR populations. Additional discussion in support of the "Dip model" can 
be found irP3 , 

In summary, the identification of the transition energy that separates galactic and 
extra-galactic cosmic rays is a central problem for high energy astrophysics, with 
wide and important implications. 

Luckily it is likely that the problem will be solved with new observations espe- 
cially those related to the mass composition and angular distribution of the cosmic 
rays. Essentially all CR models predict that the two components have very differ- 
ent mass composition (with a galactic component dominated by heavy nuclei, and 
an extragalactic one rich in protons), therefore the measurement of energy spectra 
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separately for different mass ranges lias the potential to clearly disentangle the two 
components. Similarly, the large scale angular distribution of CR below and above 
the transition energy, (possibly also for different mass groups) can allow to separate 
an isotropic extra-galactic component from a galactic one that should show the effects 
of the Milky Way magnetic confinement. 

4.4- Energy Balance for Extra-Galactic CR 

The average energy density of extra-galactic cosmic rays pl^^"^^ is the result of the 
continuous injection of particles during cosmological history, taking into account the 
losses due to energy redshift and other interactions. Integrating over all energies and 
considering only redshift losses one has: 



where Cct{z) is the power density (in comoving coordinates) of the CR sources at the 
epoch z (and / — 1 3 is an adimensional quantity). Includig pair and pion pro- 
duction lossses for protons require a more complicated but straightforward analysis. 

The estimate of the necessary CR power density is a very important constraint for 
the possible sources. This estimate is however difficult because only a small portion of 
the extragalactic spectrum is visible, with the low energy part hidden by the galactic 
foreground, and therefore the total energy density is not directly measured and can 
at best be obtained with an extrapolation based on theory. 

It is possible to estimate the energy density (and correspondingly the power den- 
sity) of cosmic rays above a minimum energy -Emin- Fig- E] plots the power density 
of the extragalactic proton sources that correspond to the fluxes given in fig. HI All 
these fluxes reproduce the highest energy part (above the "Ankle" ) of the CR spec- 
trum, but having different slopes they have different behaviours at lower energy. The 
required power density clearly depends on the minimum energy £^min, and is also a 
(rapidly varying) function of the slope a. 

The power density needed to generate the CR above 10^^ eV is of order: 



This power estimate depends only weakly on the assumed shape of the injection spec- 
trum, however extrapolating to lower energy the power density increases slowly (loga- 
rithmically) for a ~ 2 but faster (E^^) for softer injection. In the "Dip Model" where 
the CR injection is relatively soft {a > 2.4) and the minimum energy is -Emin ~ 10^^ eV, 
one has a minimum power requirement C% ^ 10^^ erg/(Mpc^s), that is significantly 
larger than the previous estimate. If the power law form injection continues with- 
out break to lower energy, the power requirement increases with decreasing -Emin and 
eventually becomes extraordinarily large. 
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Figure 6: Power density necessary to produce the extra-galactic cosmic rays, plotted as a function 
of the minimum energy -Bmin of the injected particles. 

What are the possible candidates for the sources of the extragalactic CR, based 
on this power requirement? The two leading candidates are Active Galactic Nuclei 
(AGN) and Gamma Ray Bursts (GRB). The bolometric power density of the ensemble 
of AGN is uncertain but for example it was estimated by UedaP^ multiplying by a 
factor of 30 the observed energy density in hard X-ray ([10,20] KeV) band: 

/-bolomotrix ^ O 1 n^O ( \ /r,. n 

It is not unreasonable to expect that a fraction of order 1-10% (or possibly more) 
of the released energy is transformed into relativistic particles, and therefore AGN 
satisy the energy constraint. 

Active Galactic Nuclei are the central regions of Galaxies that emit radiation in 
a very broad energy range, their source of energy is modeled as due the gravitational 
energy released in mass accretion on supermassive Black Holes (with mass as large as 
~ 10^-10^° Mq). It can be interesting to note that the power output of the ensemble 
of AGN can be checked against the observed mass density in supermassive Black 
Holes in the universe. When a mass element m falls into a Black Hole (BH) of mass 
M,, the BH mass increases by an amount AM, ~ (1 — e) m, while the energy em is 
radiated away. The relation between a BH mass and the energy radiated during its 
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formation can therefore be estimated as: 



M, — -Eradiated (22) 

The radiation efficiency e can be roughly estimated as the kinetic energy per unit mass 
{G M,/r) gained by a particle falling into the BH potential down to a radius r equal 
to a few Schwarzschild radii {Rs = 2G M,), where the particle falls is temporarily ar- 
rested by the collision with the other accreting matter. For r ^ 5 i?5 one finds e ~ 0.1. 
The measured density of BH's in the local univers^^Dl jg 

~4.6^1ixlO5M0/Mpc^ 

the corresponding density in radiated energy is of order ~ 10~^ eV/cm^, and if a 
fraction of few percent of this energy is converted into relativistic particles, this is 
sufficient to form the CR extragalactic spectrum. 

The possibility to obtain sufficient power from Gamma Ray Bursts to fuel the 
observ ed extr a-galactic cosmic ray density is the object of some dispute. Several 
authors 22]33l have argued in favor of this explanation, while according to others (see 



for exampl(p9l) the GRB energy output is insufficient. 

The long duration GRB's have been solidly associated with a subset of SN explo- 
sion. The SN blast waves are considered as the source of most of the galactic CR. 
The rate of SN explosions in the near universe is estimateci^ as: 



^observed ^ ^ Q+bA ^ (Mpc"' yr-^) (23) 

in good agreement with theoretical estimates based on the star formation (and death) 
rates. This SN rate corresponds to the power density (in the form of kinetic energy 
of the ejecta) £g" ~ 3 x 10*^° erg/(Mpc^s), and presumably to a CR power density 
5-10 times smaller. It is however assumed that most of the CR produced around 
the spherical SN shocks cannot reach the energy range where the extragalactic CR 
are visible. The GRB phenomenon, that is the emission of ultrarelativistic jets is 
associated with a small fraction of the GRB, estimated by Woosley and Howe^ as 
less than 1% the core collapse SN. Estimating a CR emission energy per GRB of order 
gj,g thjg f^Q^j^ provide a CR injection power density of order 10^^-10^^ erg/ (Mpc^s). 
This is sufficient to generate the extragalactic CR in models where the injection has a 
flat spectrum, or where the CR emission is limited to very high energy, but it is likely 
to be too small for other models with softer injection, and an emission that grows to 
E as small as a TeV or less. 



5. Cosmic Ray Point Sources 

At a sufficiently high energy the magnetic deviations of charged particles in the 
propagation from their sources should become sufficiently small so that it becomes 
possible to perform high energy cosmic ray astronomy. This simple and very attractive 
idea has been for a long time a "dream" for CR physicists and, with a high degree of 
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confidence, has recently become a reality thanks to new data of Aiiger^^'^, even if 
the result has received only weak or no confirmation by the YakutskP^ and HiReJ^ 
detectors, and the interpretation of the result remains controversial, at the center of 
very lively discussion. 




Figure 7: Hammer-AitofF (equal area) projection of the Auger (black circles) events with energy 
E > 6x 10^^ eV. The points are galaxies with z < 0.015. The thick (red) line is the so called "Super- 
Galactic plane" . The thin lines delimit regions in the sky that have received equal exposures. 

For quasi-linear propagation, the angular deviation of a charged particle propa- 
gating for a distance d in extragalactic space can be estimated assuming that the field 
has a typical value B, and it is organized in "magnetic domains" of linear size 
where the field direction roughly parallel. Summing the deviations contributed by all 
magnetic domains crossed by a particle results in the estimate of the deviation: 

^,or^N.^^....[^J^^^^<^ (24) 

(where = E/{ZeB) is the Larmor radius). Numerically: 

^20 



I /10^^eV\ 

56 = 0.53° Z vn^Mpc (— ^ j ^2^) 

Astronomy with CR of energy E is therefore possible (for an angular resolution of 
order 69) only within a sphere of radius: 

i?imagmg(^,56') = p— ^yyy- (26) 

The "CR imaging radius" is strongly dependent (oc E'^/Z'^) from the charge and 
energy of the particle considered, and shrinks also rapidly (oc (5^)^) if one requires 
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sharper images. The crucial parameter that fixes the size of the CR imaging sphere 
is the combination (B'^Xb)^^ that characterize s the e xtragalactic magnetic field. 

At the end of 2007 the Auger coUaboratiorPSESJ j^^^g published some potentially 
very important results about the angular distribution of their highest energy events. 
The Auger group has shown that the arrival directions of the 26 highest energy events 
(E'min ~ 57 X 10^^ eV), and the positions of the closest 292 Active Galactic Nuclei 
contained in the Veron Getty catalogue, and in the detector field of view with a 
maximum redshift 2;max = 0.017 (corresponding to c? < 71 Mpc) are correlated. Of 
the 26 selected events, 20 arrive within a cone of aperture ip = 3.2° around one of the 
AGN positions. The expected number of coincidences for an isotropic distribution of 
the GR arrival direction is estimated as 5.6, and the significance of the excess, taking 
into account the "statistical cost" of optimizing the three quantities {ii^min, ^^max, '?/'} 
is estimateciB as P = 1.7 x 10~^. 

This important result has attracted considerable attention an several comments 
have already appeared in the literature. The most straightforward (or "naive") inter- 
pretation of the results is of course that the AGN themselves are the sources of the 
highest energy GR, and that the angular scale of 3.2° is simply the typical deviation 
suffered by the particles during their propagation from their sources. This would be 
an extraordinary result that at the same time (a) solves (or goes a very long way 
toward solving) the problem of the UHEGR sources; (b) gives a crucial measurement 
of the properties of the extragalactic magnetic field; and more controversially (c) mea- 
sure the electric charge {Z ~ 1) of the UHEGR. It should be noted that the angular 
resolution of the Auger detector is better that 1°, and therefore the 3 degrees cone is 
not a physical, and not an instrumental effect. This deviation is accumulated during 
the entire trajectory of the particle, and can be decomposed in contributions from the 
"source envelope", the extragalactic propagation and finally deviations in the Milky 
Way disk and magnetic halo. 

{S9) ~ {69) 

source extra © {SO)uW halo © {S6)mW disk (27) 

Our knowledge of the large scale of the magnetic field of our own Galaxy especially 
away from the disk is incomplete and non-precise, however it is reasonable to expect 
that for the energies considered the integral of the magnetic field / di B{i) over most 
line of sights corresponds to deviations larger than 3° if the charge much larger than 
unity, and therefore one can use our own galaxy as a spectrometer to determine that 
the GR are protons or at most light nuclei. Interpreting the dimension of the cone 
that optimizes the correlation with the AGN as the deviation angle of the particles 
one can also use equation (IMl) . for a first estimate of B'^ Xb- Inserting the values 

'^This probability is estimated in robust way, "optimizing" the three parameters from only the 
first half of the collected data, and then computing P from the second half of the data, keeping the 
parameters fixed 
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59 = 3.r, ~ 71 and E = 60 X 10^^ eV in one obtains: 

5^ \b - 0.15 (nGauss)^^ Mpc (28) 
(where Z ~ 1 is the charge of the observed particles). This imphes an imaging radius: 

= (l/^)' MP-^ (29) 

that is very encouraging for the perspectives of CR astronomy. 

These "nominal" interpretation of the Auger data has however a number of signif- 
icant problem and it can only considered as tentative. In fact the publications of the 
Auger group are very careful in observe that since the AGN distribution is correlated 
to the dis tributi on of normal galaxies. Moreover, the Auger data can can also be 
interpretecPSHOJ assuming that the anisotropy the effect is mostly due to the a the 
contribution of a single source, the nearby {d ~ 3.5 Mpc) Active Galactic Nucleus 
GEN A, generating as much as one third of the events above 6 x 10^^ eV, with a 
larger angular spread of order 69 ~ 10° (or more). Substituting the shorter distance 
and the larger 69 one obtains an estimate of Xb larger by a factor of order 200: 

B"^ \b ~ 42 (nGauss)^ Mpc (30) 

and correspondingly a 200 times smaller imaging radius. The discrimination between 
the two estimates (l28l) and (1301) for the parameter that describes the extragalactic 
magnetic field is obviously crucial for the future of the field. 

6. Cosmic Ray Composition and Hadronic Interactions 

A crucial problem for GR science is the estimate of the energy and mass of the 
detected particles. At high energy it is only possible to observe the showers produced 
in the atmosphere by GR particles, and to extract information about the properties of 
the primary particle it is obviously necessary to have a sufficiently accurate description 
of the shower development. The key ingredient, and the dominant source of systematic 
uncertainty is the description of the properties of hadronic interactions, such as: cross 
sections, final state multiplicities and inclusive energy spectra. At a fundamental 
level there are few doubts that Quantum-Ghromo-Dynamics (QGD) is a complete 
and successful theory of the strong interactions; however the fundamental fields that 
enter the QGD Lagrangian are only quarks and gluons, and perturbative calculations 
are only possible for high momentum transfer processes between these fields. Most 
of the observables that are relevant for shower development cannot be calculated 
from first principles. These quantities (like cross sections and multiplicities) must 
be obtained from measurements obtained in accelerator experiments. The difficulty 
here is that the available data do not entirely cover the relevant phase space, and 
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Figure 8: The left panel shows the angular distribution of the 26 highest energy events of AugeJ^^. 
The horizontal x axis corresponds to the right ascension and the vertical x axis to the declination 
i5. The y coordinates have been "stretched" taking into account the declination dependence of the 
exposure, so that for an isotropic distribution, the points that correspond to the measured directions 
should fill uniformly the square. The probability that the measured arrival directions are distributed 
uniformly is few percent. The right panel shows the distribution of near {z < 0.015) AGN in the 
Veron-Cetty catalogue. The position of CEN A is also shown as a red circle in the left panel. 

therefore for the interpretatfon of CR observations it is necessary to extrapolate from 
the regions were observations are available, using guidance from the known structure 
of the underlying theory. 

The incomplete phase space coverage is evident if we consider the center of mass 
energy of the first CR interaction, that for the highest energy particles with E ~ 
10^° eV is a/s ~ 430 TeV, four hundred times larger than the highest energy where 
accelerator data is available (at the pp Tevatron collider). This "energy gap" will be 
only partially filled by the eagerly awaited data of the LHC collider at CERN with 
a cm. energy ^/s = 14 TeV (see for example a plot of the total pp cross section in 
fig. [9]). It should also be noted that data at high a/s is available only for pp or pp 
collisions, (and not for meson-nucleon and hadron-nucleus interactions). Note also 
that for an accurate description of the shower development it is necessary to have a 
precise description of the highest energy secondaries (in the projectile fragmentation 
region), these particles are experimentally the most difficult to measure in collider 
experiments, and significant uncertainties about their energy spectrum still exist. For 
this reason even below the knee {^/s ~ 2.3 TeV) uncertaintis in the description of 
the hadronic interactions are a source of significant errors in the energy and mass 
measurement. 

Fortunately, the uncertainties in the description of hadronic interactions play only 
a limited role in the energy determination for the highest energy particles, where one 
can observe showers using the fiuorescence light method pioneered by the Fly's Eye 
detector. Relativistic charged particles in air excite nitrogen molecules in the medium. 
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Figure 9: Total cross sections for pp and pp scattering. The filled (empty) points are measurements 
for the pp (pp) channel. The solid line is a is the result of a fit in ^Jj The dotted line is the fit in 

m 



These excited molecules can return to the fundamental states emitting isotropically 
fluorescence photons, that are detectable at the ground, with the shower appearing 
as a quasi point-source of light moving at the speed of light across the sky. The 
photon emission is proportional to the number of charged particles, and thereore 
from the angular distribution of the detected photons it is possible to reconstruct the 
longitudinal profile Ne{X) of each detected showeiQ. The longitudinal profile Nf,{X) 
allows the model independent reconstruction of the energy dissipated in air by the 
shower as ionization (and then ultimately heat): 

/■^ground / dE \ 

^ionization = I dXl—)N,{X) (31) 

The total shower energy can be calculated adding to this dominant contribution 
smaller correction terms that take into account the energy that goes into neutrinos, 
muons, and is dissipated in other (as hadronic and electromagnetic components) below 
the ground: 

-^shower -^ionization ~l~ Ey -\- -\- -Egj-ound ('^^) 

' Ne is the number of charged particles present at the depth X, and the column density X can 
be used as a coordinate along the shower axis 
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The correction terms E,^, and -Egmund are mass and model dependent, but fortu- 
nately account for a small fraction (of order ~ 10%) of the primary particle energy, 
and the energy measurement has only a small model dependence. 

In other words, according to equation ( l3Ti) . the area subtended by the longitudinal 
profile iVe(X) is a measurement of the shower energy that is independent from the 
particle type and from the modeling of hadronic interactions; on the other hand the 
shape of the curve depends on both the primary particle mass and the properties of 
hadronic interactions. 

The shape of the showers longitudinal development is in fact the best available 
method to estimate the mass of UHECR. It is in fact intuitive that, for the same 
total energy, the showers generated by a nucleus of mass A will develop more rapidly, 
reaching their maximum size at a shallower depth. Because of fluctutations this 
method cannot identify the mass of an individual shower, but can be used statistically 
to estimate the mass composition at a given energy. The simplest method is to use 
the average value (Xmax (-£')) of the position of maximum size for all showers with 
energy ii^ as a measurement of the average mass of CR at that energy. 

For a more quantitative analysis one can denote X'^^J^E) the average position of 
shower maximum for proton showers and, since showers penetration grows approx- 
imately logarithmically with increasing energy, develop in first order around energy 
E*. Measuring the energy in units of E* one obtains: 

Xl,^{E)^X,{E*)+D, InE (33) 

The quantity Dp gives the increase in the average position of shower maximum when 
the p energy increases by a factor e and is known as the proton "elongation rate" 

To a very good approximation the quantity X^^^{E), that is the average position 
of maximum for showers generated by a primary nucleus of mass A and total energy 
E is related to the average maximum for proton showers by the relation: 

XtAE) ^ X^ax (I) ^ X^ax(^) -Dp\nA (34) 

In ( !34l) the first equation is clearly correct for a naive superposition model (where the 
shower generated by a nucleus containing A nucleons is described as the superposition 
of A nucleon showers each having energy E /A), but in fact it has a much more general 
validitjJ^, since it is based on the relations between the cross sections for nuclei and 
nucleons in the framework of the Glauber model^; the second equality in flM|) follows 
from (1331). 

More in general, if the CR of energy E are a mixture of nuclei with different mass 
A, the average (Xmax) can be expressed as: 

{X^UE)) ^ Xl,^{E) - Dp{E) (In A)^ ^X. + D, [\^E - (In A)^] (35) 



27 



where (InA)^ the average value the logarithm of A for particles of energy E. One 
can then use equation fl35l) to estimate the mass composition of UHECR comparing 
the measured (Xmax) with the value expected theoretically for protons: 



{lnA)E 



Dp 



(36) 



Similarly comparing the experimental elongation rate -Dexp = d{X^i^^{E)) / dlnE with 
the elongation rate for protons (or nuclei) Dp, one can estimate the rate of variation 
of the composition with energy: 



An illustration of the present situation of the mass composition of UHECR is 
shown in fig. [TDl The points in the fi gur e show the average X^ax for UHE showers 
obtained by the HiRes-MIASSl, HiRes^ and AugeJ^ experiments. The thick lines 
are fit to the data points obtained by the experimental groups, while the thin lines in 
the figure are theoretical predictions calculated assuming a pure composition of only 
protons or only iron nuclei based on the hadronic models of SibylP^ and QGS JEtS3F 
At face value these results seems to indicate that the UHECR are not pure protons 
nor pure iron nuclei, but are a combination (with approximately equal weight) of 
the two species, or are mostly nuclei of intermediate mass. This conclusion however 
relies on the assumptions that the models used for the description of the hadronic 
interaction are at least approximately correct. 

The composition of the CR has very important consequences for the properties 
of their acceleration sites, and therefore it is of great importance to gain a more 
robust understanding of the relevant properties of hadronic interactions, and a better 
estimate of the systematic uncertainties. It is clear for example that if the cross 
section is larger or if the energy spectra of the final state particles are softer than 
what is assumed in the theoretical prediction, the real showers develop faster than 
the simulated ones, and the composition interpretation is biased in the direction of a 
too heavy composition. 

In the discussion of the previous paragraph, information about elementary par- 
ticle physics (obtained from accelerator experiments) is used to measure the mass 
composition of the highest energy CR, and then infer the structure and properties of 
their astrophysical sources. It is interesting that this "logical flow" can be in principle 
be "inverted" . In fact there are (at least in principle) methods to measure the com- 
position of CR that are independent from the detailed properties of their showers. 
These methods, as discussed above, are based on essentially two ideas: (i) the "cosmic 
magnetic spectrometer", and (ii) the imprints of energy losses on a smooth injection 
sprectrum. If the "nominal interpretation" of the Auger anisotropy analysiJ^ is 
correct and the deviations of extragalactic particles with ~ 6 x 10^^ eV is as small 



d{\iaA)E 
d\nE 



1 - 




exp 



(37) 
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Figure 10: Measurements of (X^ax) obtained by the HiRes-MIA HiRes SH and Auge^ 
experiments. The lines are monteca rlo c alculations for protons (jj) and iron nuclei (Fe) performed 
using the SIBYLL '431 and QGSJEt'44^ interaction models. 



as 3° their electric charge is in fact strongly costrained, by our (however still poor) 
knowledge of the galactic magnetic fields. Similarly, if the "Dip Model" of Berezin- 
sky and coUaboratorJ^Sl jg correct and the structure of the ankle corresponds to the 
effects of e~^e~ pairs production in py interactions, then the identity of the UHECR 
particles is determined. In this case, the requirement of consistency for the pene- 
tration of the corresponding showers can be used to obtain information about the 
properties of hadronic interactions. 



7. Outlook 



The observation and understanding of the "High Energy Universe" is one of the 
important frontiers of fundamental Science. Our universe contains astrophysical ob- 
jects that are capable to accelerate particles to ultrarelativistic energies as large as 
~ 10^° eV. The nature and structure of these cosmic accelerators and the physical 
mechanism that are operating in them are beginning to be clarified, even if large uncer- 
tainties still exist. In fact, several classes of these cosmic accelerators are now known 
to exist. There is good evidence that hadronic particles (protons and nuclei) are accel- 
erated in SuperNova Remnants and Active Galactic Nuclei, while electrons/positrons 
are also efficiently accelerated in several other environments like pulsars and pulsars 
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wind nebulae. Gamma Ray Bursts are also clearly sites of particle acceleration. These 
extraordinary environments are very important astrophysical laboratories to test our 
physical theories. 

The acceleration of charged particles is unavoidably connected to the radiation of 
photons and (in the case of hadronic particles) to the emission of neutrinos. More- 
over many of these accelerators are connected to the violent motion of large masses 
and to the possible emission of gravitational waves. Therefore the complete study of 
these objects must rely on the combined use of different "messengers" (photons in 
a very broad range of wavelength, cosmic rays, neutrinos and gravitational waves). 
We are now witnessing the opening of new "astronomies" with cosmic ray parti- 
cles (at sufficiently high energy to reduce the bending by magnetic fields), neutrinos 
and gravitational waves. The nearby active galactic nucleus Cen A is very likely 
the first astrophysical objects imaged with charged particles. The development of 
multi-messenger astronomy is a revolution comparable to the use of the telescope for 
astronomical observations started by Galileo in 1609, nearly exactly four centuries 
ago, and is certainly going to produce exciting results. 
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